home *** CD-ROM | disk | FTP | other *** search
/ AmigActive 24 / AACD 24.iso / AACD / Graphics / vcdimager / source / differ2.txt < prev    next >
Text File  |  2001-07-02  |  3KB  |  105 lines
  1. diff -Pc3prb ram:old/rint.c ./rint.c
  2. *** ram:old/rint.c    Thu Jan  1 00:00:00 1970
  3. --- ./rint.c    Mon Jul  2 19:15:43 2001
  4. ***************
  5. *** 0 ****
  6. --- 1,98 ----
  7. + 
  8. + 
  9. + #include <math.h>
  10. + 
  11. + #if defined(vax)||defined(tahoe)
  12. + 
  13. + // Deal with different ways to concatenate in cpp
  14. + #  ifdef __STDC__
  15. + #    define    cat3(a,b,c) a ## b ## c
  16. + #  else
  17. + #    define    cat3(a,b,c) a/**/b/**/c
  18. + #  endif
  19. + 
  20. + // Deal with vax/tahoe byte order issues
  21. + #  ifdef vax
  22. + #    define    cat3t(a,b,c) cat3(a,b,c)
  23. + #  else
  24. + #    define    cat3t(a,b,c) cat3(a,c,b)
  25. + #  endif
  26. + 
  27. + #  define vccast(name) (*(const double *)(cat3(name,,x)))
  28. + 
  29. + 
  30. + #  define vc(name, value, x1,x2,x3,x4, bexp, xval) \
  31. +     const static long cat3(name,,x)[] = {cat3t(0x,x1,x2), cat3t(0x,x3,x4)};
  32. + 
  33. + #  define ic(name, value, bexp, xval) ;
  34. + 
  35. + #else    //  vax or tahoe
  36. + 
  37. +    // Hooray, we have an IEEE machine
  38. + #  undef vccast
  39. + #  define vc(name, value, x1,x2,x3,x4, bexp, xval) ;
  40. + 
  41. + #  define ic(name, value, bexp, xval) \
  42. +     const static double name = value;
  43. + 
  44. + #endif    /* defined(vax)||defined(tahoe) */
  45. + 
  46. + 
  47. + 
  48. + vc(L, 4503599627370496.0E0 ,0000,5c00,0000,0000, 55, 1.0) /* 2**55 */
  49. + 
  50. + ic(L, 4503599627370496.0E0, 52, 1.0)              /* 2**52 */
  51. + 
  52. + #ifdef vccast
  53. + #define    L    vccast(L)
  54. + #endif
  55. + 
  56. + 
  57. + #if defined(vax)||defined(tahoe)      /* VAX D format */
  58. + #include <errno.h>
  59. +     static const unsigned short msign=0x7fff , mexp =0x7f80 ;
  60. +     static const short  prep1=57, gap=7, bias=129           ;
  61. +     static const double novf=1.7E38, nunf=3.0E-39, zero=0.0 ;
  62. + #else    /* defined(vax)||defined(tahoe) */
  63. +     static const unsigned short msign=0x7fff, mexp =0x7ff0  ;
  64. +     static const short prep1=54, gap=4, bias=1023           ;
  65. +     static const double novf=1.7E308, nunf=3.0E-308,zero=0.0;
  66. + #endif    /* defined(vax)||defined(tahoe) */
  67. + 
  68. + 
  69. + double copysign(double x,double y)
  70. + {
  71. + #ifdef national
  72. +         unsigned short  *px=(unsigned short *) &x+3,
  73. +                         *py=(unsigned short *) &y+3;
  74. + #else    /* national */
  75. +         unsigned short  *px=(unsigned short *) &x,
  76. +                         *py=(unsigned short *) &y;
  77. + #endif    /* national */
  78. + 
  79. + #if defined(vax)||defined(tahoe)
  80. +         if ( (*px & mexp) == 0 ) return(x);
  81. + #endif    /* defined(vax)||defined(tahoe) */
  82. + 
  83. +         *px = ( *px & msign ) | ( *py & ~msign );
  84. +         return(x);
  85. + }
  86. + 
  87. + 
  88. + 
  89. + double rint(double x)
  90. + {
  91. +     double s,t;
  92. +     const double one = 1.0;
  93. + 
  94. + 
  95. +     if (x != x)                /* NaN */
  96. +         return (x);
  97. + 
  98. +     if (copysign(x,one) >= L)        /* already an integer */
  99. +         return (x);
  100. +     s = copysign(L,x);
  101. +     t = x + s;                /* x+s rounded to integer */
  102. +     return (t - s);
  103. + }
  104. + 
  105.